Manuscript for Revista Mexicana de Astronomia y Astrofisica (2007) 



THE EFFECT OF SUBLIMATION TEMPERATURE 
DEPENDENCIES ON DISK WALLS AROUND T 
TAURI STARS 

Erick Nagel/ Paola D'Alessio,^ Nuria Calvet,^ Catherine Espaillat,^ and 
Miguel Angel Trinidad^ 
Draft version: October 9, 2012 

RESUMEN 

El polvo no puede sobrevivir mas cerca de la estrella de un punto donde al- 
canza una temperatura igual a la temperatura de sublimacion. La frontcra 
entre una region sin-polvo y polvosa define la pared de sublimacion. En la lit- 
eratura dos modelos para la estructura de la pared son usados: una pared con 
una temperatura de sublimacion fija y con una temperatura de sublimacion 
^ dependiente de densidad. En la primera, la pared es vertical y en la segunda, 

^ J la pared es curva. Encontramos diferencias importantes entre la SEDs de es- 

tos modelos en el intervalo de longitudes dc onda desde 3 a 8/xm, sicndo la 
emision de la primera mas grande que la segunda. Cuantificamos las diferen- 
cias en los colorcs dc IRAC cuando cstos modelos son usados para explicar las 
observaciones. En el diagrama de IRAC color-color ([3.6]-[4.5] vs. [5.8]-[8.0]), 
los modelos cstan localizados en rcgioncs cspecfficas, dada la inclinacion, la 
tasa de acrecion de masa, o cual de los modelos es usado. 

ABSTRACT 

The dust cannot survive closer to the star from the point where a grain reaches 
a temperature equal to the sublimation temperature. The boundary between a 
dust-free and a dusty region defines the sublimation wall. In the literature two 
models for the structure of the wall are used: a wall with a fixed sublimation 
temperature and a wall with a density-dependent sublimation temperature. 
In the former, the wall is vertical and in the latter, the wall is curved. We 
find important differences between these models SEDs in the wavelength range 
from 3 to 8/^m, being the emission of the former larger than that of the latter 
model. We quantify the differences in IRAC colors when these models are 
used to explain the observations. In the IRAC color-color diagram ([3.6]-[4.5] 
vs. [5.8]-[8.0]), the models are located in specific regions, either depending on 
the inclination, the mass accretion rate, or which model is used. 
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1. INTRODUCTION 
One of the motivations for the study of the dust in the proto planetary disks 



inner regions is the fact tha t terrestrial planets are formed there (jGoldreich fc Ward 



197i lAlibert et al~ll2010l) . The dust might not survive if the dust grains 



reaches temperatures higher than the sublimation temperature. Thus, the 
innermost region of a disk around a star is dust free and has an opacity hole 
or deficit. The outer boundary of the hole formed in this way is called a sub- 
limation wall or the disk inner rim. This surface separates an outer dusty and 
an inner gaseous disks, and its shape depends on the characteristics of the gas 
and the dust. 



Hillenbrand et al. I (|1992I) model the near-infrared (NIR) excess of Hcrbig 
Ae/Be stars as coming from material located at the inner rim of the disk with 
a temperature around 1500K. This is consistent with the evaporation tem- 
perature of silicate grains, thus, they interpret that the emission is produced 
in a sublimation wall. In fact, the edge of the hole produced by dust sublima- 
tion mostly emits in the NIR. The modelin g of the NIR excesses as produced 



in this kind of walls was followed by many ( Tuthill. Monnier fc Danchi 2001 : 

Eisner et al. boOStlAkeson et al. boOStlMonnier et al. IboOStlMihan-Gabet et"aL 



20071 ). A nice characterization of the shape of the wall using a Monte Carlo ra- 



dia tive transfer cod e , takin g into account stratification of grains was developed 



by IWhitnev et al. I (|2004l ) . They found an empirical formula for the radius 
of the sublimation wall {Rwaii) in terms of the evaporation temperature Tsub 
{Rwaii oc which is consistent with simplified expressions commonly 

used. The range of dust sublimation temperature in terms of gas density 
is between 10~^^ an d lO^^.gcm^'^ and for different grain species is given by 



Pollack et al. I (|1994[ ). For the sihcate olivine, the range is (929 - 1774)i^; for 
the silicate pyroxene, Tsub = (920 — 1621)ii'; for troilite, Tgub = 680/^; and 
for iron the range is (835 — 1908)X. 



For the case of very low luminosity systems as the brown dwarfs . iMavne fc Harries 



(|2010l ) notes that the location of the dust sublimation wall is equal to the 



disk co-rotation rad ius with the magnetosphere, as it is also suggested by 



Eisner et al. I (|2007[ ) . For this case and also for configurations where the sub- 
limation wall is inside the co-rotation radius, the magnetic field truncates the 
disk, thus, the resulting shape of the wall depends strongly on this processes 
and weakly on the sublimation phenomenon. In our case, the stellar lumi- 
nosity is large enough to move the wall outside the magnetosphere, thus, its 
geometry is given by the physical state of dust and gas. 

Taking a unique sublimation temperature, leads to a unique radius, thus, 
a vertical wall. We refer to this model as T^.fix wall. For this model, we use 
a sublimation temperature characteristic of the disk midplane, and assume it 
is the same at every height (in spite of the variation of density with height). 
This implies that the wall is vertical, i.e., the inner surface of a hollow cylinder 
, but also, that the surface temperature is constant (and equal to the assumed 
sublimation temperature). On the other hand, taking into account the de- 
pendence of the sublimation temperature with density, the wall is curved. We 
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refer to this model as Tq.p wall. The variation of density with height and 
radius in disks models are taken from structures assumed to be in vertical 
hydrostatic equilibrium. Again, there is a twofold effect, affecting the solid 
angle of the visible portions of the wall and also, the surface temperature of 
the wall, which is different at each pixel. When we compare the SED of a 
Tojix wall with that of a To^p wall, both effects are present, i.e., differences in 
area and in temperature, but them are difficult to disentangle. 

If the gas density and the dust composition do not depend on the verti- 
cal coordinate, then a Tq ji^ wall is formed. We want to point out the fact 
that sometimes due to the unknowns of the composition and density profile, 
one is naturally leaded to simplify the system and assume a homogeneous 
vertical distribution of matter. However, as we will see in the following, a 
vertical stratified model, curves the wall . Previously, the assumption was 



that the shape of the wall was vertical (jPuUemond. Dominik fc Natta 112001 



D'Alessio et al. 1120051 ). However, for the mo deling of the stationary state of 



the disk around EX Lup ( Sipos et al. 20091) . a Tojix wall was unable to ex- 



plain the IR observations, thus an ad hoc rounded inner r i m was required. 
A non- vertical wall was also considered by llsella fc Natta I (120051 ) , they in- 
terpreted a 2/im emission bump, in observed SEDs of Herbig Ae stars, as 
coming from a sublimation wall. They consider that due to the density de- 
pendence on the sub li matio n temperature, the wall is curved. Two years later 



Tannirkulam et al. 



(120071 ) pursues this further, taking into account grain 
sedimentation. The dust is composed of two grain size distributions, char- 
acterized by different scale heights. In these w orks a Tn.^ w all model can 



reproduce the Herb i g Ae stars NIR spectrum, llsella &: Natta I (|2005l ) and 



Tannirkulam et al. I (|2007() study the differences between synthetic images of 



Tojix and To^p wall models, here, we analyse the effects on the colors. Un- 
fortunately the spatial resolution and, limited sensitivity and the number of 
telescopes of NIR interferome tric observations, make difficu lt to confirm the 



geometry of the inner region (jPullemond fc Monnier 1120101 ). Thus, we have 



some intrinsic degeneracy on the models used to interpret the data. Due to 
this, it is important to include all the physics that we can on the models of 
these walls. 

In this work, we consider that sublimation is the mechanism responsible 
to produce the inner hole, the physi cs involved in this pr ocess shapes the wall. 



Note that a recently formed planet (lOuillen et al. 1120041) and a photodissocia- 



tion fiux (IClarke. Gendrin fc Sotomavor 1120011 ) are also able to create a hole. 
The magnetic rotational instab ility is responsible to create a w ind, which 



in turn is able to form a hole (jSusuki. Muto fc Inutsuka I l2010l ). A binary 



syste m is another way to create a hole, in this case by gravitational in terac 



tions ( Hartmann et al. l2005at l lEspaillat et al. l2007l : iNagel et al. l2010l) . The 
physical peculiarities of each process will define the structure for the wall. 



In a sense, this work follows the steps of llsella fc Natta I (|2005() . because 
we are taking into account the same physics to describe the grain sublimation. 
However, one difference is that we focus on T Tauri instead of Herbig Ae stars. 
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Our first aim is to compare models of Tg^p and T^jix walls, characterizing pa- 
rameters like location and surface temperature. Another questions to address 
are: if we change the inclination or the mass accretion rate, still is possible to 
discriminate between a T^jix and a T^^p wall based only on the SED?, based 
on the Infrared Array Camera (IRAC, on board the Spitzcr Space Telescope) 
observations, does a model of a disk plus a T^jix wa ll or a disk plus a T^ n 
wall show differences in the IRAC colors? Noteworthy. Ilsella fc Natta I ( 20051 ) 
mention that based on an image is easy to discriminate between them, because 
for the Tojix one, only the half part of the wall that is farthest away from the 
observer show emission in the line of sight. The observed To^p wall emission 
comes from every azimuthal angle of the wall. Without an image, one cannot 
choose between both models, however, looking the differences in the SED is a 
way to argue if for the problem at hand it is enough to take a Tojix model. 

The minimum location of the wall (where the maximum temperature oc- 
curs) is a parameter that changes when using Tojix or Tg^p wall models, one 
can conclude in the following sections that this parameter differ at most by 
10% (see Table[T]). In the case of the IRAC colors, between the Tg^p and T^jix 
wall models, the [3.6]-[4.5] and [5.8]-[8.0] colors changes in around 20% and 
10%, respectively. In §[2] we present the details of the code used in this work, 
followed in §[3] by the resulting characteristics of the walls, cither the Tojix 
l3.ip or the To^p walls 13. 2p . §|4] shows the IRAC colors associated to the 
models presented. Finally, in §[5]wc present a summary and the conclusions. 



2. DESCRIPTION OF THE CODE 



The To fix wall emission is calculated with the code used in D'Alessio et al 



(|2005[ ) for an isolated star or for a binary system in iNagel et al. I (|2010h . We 



consider that the wall is optically thick but the emission from an optically 
thin atmosphere is taken into account. The wall is heated by the impinging 
geometrically diluted stellar radiation flux coming from the photosphere of 
the star and from the accretion shocks. We assume that the stellar radiation 
is plane parallel. Thus, for a Tojix wall, the radiation arriving at each point is 
the same. The total emission is the addition of the contribution of each layer 
at given r, extingui shed with the niateri al in front of it. The temper ature is 



calculated following ICalvet et al. I (|l99lh and iD'Alessio et al. I (|2005l ) . 



We have a. s sumed that the opac i ties ar e ind ependent of t, in a ccord ance of 



Calvet et al. I (|l99lh)Calvet et al. I (|l992l ) and ID'Alessio et al. I (|2005[ ). This 



is a necessary assumption in order to find an analytical expression for T{t). 
The temperature at every depth of the wall atmosphere is lower than the 
sublimation temperature. Thus, there is no sublimation of dust, and it is safe 
to take a constant opacity. 

The geometrical effects producing shadowing of regions of the wall by 
regions of the wall closer to the observer are taken into account. The SED 
of the wall is calculated integrating the flux of every point in the wall whose 
normal has a component in the direction of the observer, times the solid angle 
subtended by every pixel. 
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The minimum {a„iin) and mciximum size (ctmaa:) of the grains are 0.005^m 
and 0.25/xm, respectively; the power law exponent is —3 .5; which are param- 



eters typical for interstellar grains (jMathis et al. 1 119771 ). The dust is com 



posed of silicates (pyroxenes, Mgo.s Feo.2 SiO^; and olivines, Mg Fe SiO^ ), 
graphite and troilite. We adopt a dust-to-gas mass ratio for the silicates, 
Csi! = 0.0034 lOraine fc Lee 1119841 ): for the graphite, Cgrap = 0.0025; and for 
the troilite, C,troi = 7.68 x 10~^. The com position and abundance are typical 
for accretion dis ks ( Pollack et al. Il994) • The optical properties of the sili- 



cates come from iDorschner et al. I (|l995l ). From the dust composition chosen, 
the silicates are the grains with the highest sublimation temperature. From 
this follows the fact that the silicates rules the location and shape of the wall. 
Dust species with a higher sublimation temperature as corundum (74^203), 
in principle will affect the location and stru cture of the wall, beca use such 



grains are the ones formed closest to the star (jVerhoelst et al. 1120061 ). A con- 



sistent modeling of a disk with corundum requires a study that at the same 
time takes into account the gas and dust opacity, because in the temperature 
range where corundum is formed, the gas contribution to the opacity is a siz- 
able fraction of the total opacity. In the region where the silicates grains are 
formed, the silicates opacity is arou nd 6 orders of magnitude larger than the 



gas opacity (jFerguson et al. 1 120051 ) . thus, it is not necessary to include the 



gas contribution to the case treated here. This is the reason why in the sub- 
limation wall formation is not important to know the gas opacity, something 
that we cannot leave aside when corundum is included in the mixture. 

A ncar-IR emission study of a disk-star system should include the con- 
tribution of a gaseous disk inside the sublimation wall, the importance of 
this is highlighted i n the inter pretation of interfero metric observations by 



Tannirkulam et al. 



(|2008l ) and Eisner et al. I (|2010l ). Either for the study 



of the dust-free gaseous disk emission or the shaping of the wall by the pres- 
ence of corundum, a detailed knowledge of the gas opacity is necessary. This 
is a non trivial issue, because of the presence of millions of lines and also, 
because it is not clear what kind of mean opacity is representative for the ap- 
proach used here. Due to this, the gas emission problem is beyond the scope 
of this paper, but should be taken into account in the future. 

The Tp.p wall emission is calculated based on the code just described, but 
including an arbitrary shape for the wall. Unlike a T^jix wall, for the Tg^p 
wall, the radiation arriving to different places of the wall is not the same. 
The impinging flux depends on the angle of incidence a, (the angle between 
the normal to the wall surface and the incidence ray), specifically is propor- 
tional to cos a, which in turn depends on the wall shape. Thus, we have to 
characterize a in order to get a value for the temperature in the wall. In 
other words, we need to know beforehand a to get the temperature, but we 
require the temperature to calculate a, this means the wall shape. A way to 
solve this problem is to note that if the scattering of the stellar radiation is 
neglected along with the heating from inner regions (viscously produced), an 
expression for T(t = 0) without dependence on a is found. The contribution 
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of the scattered emission has the characteristic frequency range of the stellar 
radiation, thus, its main contribution is at wavelengths smaller than the peak 
of the emission of sublimation walls. Because of this, to neglect the scatter- 
ing is reasonable when one focus in the infrared frequency range. Using this 
expression, the wall shape is defined as the points where this temperature 
is equal to the sublimation temperature. This turns out in a shape for the 
wall surface; from where an incidence angle can be calculated. With these 
assumptions, the temperature as a function of r can be written as: 

rpf ^4 {L^+ Lace) COS a , 

T{t) = — ci+C2e --°), (1) 



where q = Kinc/nd, 



and 



3 cos a 

ci = , (2) 



q 3 COS a 
C2 = ■ (3) 

COS a q 



Then, the temperature at the surface (r — 0) is 



4 (i^ + Lace) cos a / 3 cos a / q 3 cos a 



lotJTrr^ \ q ycosa 1 / J 

and we k now that it should be equal to the sublimation temperature. 



Tsub{p{z, R)) (jPollack et al. Ill994[ ). which depends on density. This results in 
an equation for a{z,R). 

Here, Kmc and Kd are the mean opacity of true absorption, evaluated 
at the temperature characteristic of the incident stellar radiation, and at the 
temperature of the disk, respectively. L^, is the luminosity of the star and Lace 
is the luminosity produced by the shocks of the material accreting along the 
magnetic field lines. Note that q depends on the temperature, thus equation^] 
is solved itcratively. Knowing the temperature, the radius of dust destruction 
is calculated substituting t = in equation |4l thus. 



/ \2 ^aee ( ^ine\ 1 /-\ 

rdes^z) = —- -— j-, (5) 



16na \ Kd J Tsub{z) 
in which the scattering and the local radiation field arc neglected. 

3. CHARACTERISTICS OF WALLS 

The diflFerences in the SED when comparing a Tp^p and a T^jix wall, are 
the result of a combination of, at least, 2 effects: 1) geometry, because in a 
To.p wall, each pixel shows a different effective area to the observer than in a 
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Fig. 1. The shape of the Tojix waUs are shown with thin lines. The To,p walls shapes 
are presented with thick lines. Models presented are M = 1.625xlO~*MQyr-i (solid 
line), M = 3.25 x lO^^Mgyr"^ (pointed line), and M = 6.5 x IQ-^ M^yr-^ (dashed 
line). 



Tojix wall and 2) surface temperature gradient, because what we are assuming 
that curves the wall is the dependence of the sublimation temperature with 
density, thus this implies that at each height, the surface of the wall would have 
a different temperature than a T^ j^^ wall, defined with a unique sublimation 
temperature. 

In this section, we calculate the wall emission using the code described in 
§ [2j Our intention is to construct models for the wall emission for a typical 
young low mass star: = 0.5 M,^, = 2.0E(t,,T^ = AQQQK, and M = 



3.25 X IQ-^MQyr-^ (jGullbring et al. Ml998l ). This is our fiducial system. 



Two kinds of models are considered, the first one is a Tojix wall with constant 
surface temperature and the other is a wall, with a shape given by how the 
sublimation temperature depends on density (see Figure [1]). 

As noted in § [1] the modeling of interferometric observations primarily 
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depends on two basic parameters: a typical distance to the emission region, 
and a typical temperature. In order to compare these parameters with other 
works, we present in Table [T] for all the models, either the Tgjix or the Tg^p 
walls; the values of the minimum radius of the wall and the temperature at 
this location. For comparison, i n Table [J the parameters for the models of 
the 4 T Tauri stars presented in Eisner et al. ( 20091 ) are shown. We are not 
looking a concordance between the values. As a matter of fact, we expect 
quite the opposite because this set of values comes from two different ap- 
proaches. Our modeling produces these v alues from a detaile d description of 
the wall formation. From the other hand, lEisner et al. I (120091 ) interpret their 
interferometric observations with R y,„u a nd Tu}aii as two free parameters of a 



simplistic model. In lEisner et al. I (I2OO9I) ' model, these parameters are fitted 
without a model of the star or dust composition. 

An important fact to notice here is that to obtain the interferometric 
image of the infrared observational data is required a prior knowledge of the 
observed structure. Then the physical parameters extracted from infrared 
interferometric observations are calculated using a particular model; in other 
words, they are model dependent. These difficulties are not present in the 
radio interferometry, where the large number of points in the visibility plane, 
along with the high angular resolution and sensitivity obtained, allows to 
apply the inverse Fourier transformation to get a real image of the object, that 
does not depend on a particular model. In the case of NIR interferometry, 
the number of visibilities is so low that we cannot invert the problem, thus, 
we depend on a model to extract the parameters of the system. 

Summarizing, we should be careful to compare the set of parameters ex- 
tracted in this way w ith values given by a model independent of observations. 



Eisner et al. I (|2004l ) uses 5 different models to interpret the interferometric 
observations at 2.2^m of 14 Herbig Ae/Be stars. These span geometries such 
as an envelope, disk or ring. Looking for the best fit allows to choose the 
model, however, this do not completely guarantee that this is the right model. 
A qualitative comparison between the models presented here and the model 
consistent with observations can be done (see Table [Ij. 

For a Tb,p sublimation wall mod el, we note that the sublimation tem- 



perature {Tsub) depends on density (jPollack et al. Ill994l ). Due to the fact 



that the density decreases with height and that Taub increases with density, 
the modeled shape is convex. The lower denser parts of the wall are closer 
to the star; at high altitudes the density decreases, thus Tgub decreases and 
the wall moves further out. The density structure is given by a disk mod- 
eled with a detailed 2D numerical solution of the radiation transfer equations 

Models built with this assumption arc previously 



(D'Alessio et al. 



19981) 



given bv llsella fc Natta I (j2005l ) for Herbig Ae stars. 



3.1. T, 



OJix 



WALL 



A- Tojix wall is defined with a constant sublimation temperature in its 
surface; here we take Tsub = UOOK. For the luminosity of the typical low 
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TABLE 1 



TABLE OF MAXIMUM Twall AND MINIMUM Rwall- 



Model 


min{Rioaii){AU) 


max{Twaii){K) 


Tn f (M — 1 fi9^ X ^C\~^Mr^.llr~^\ 


n 071 5 


1 4on 


-'0,/i2;V" — O.Zd X iu IVlQi/l ) 


OTQ 

U .VJ ( fy 




ToJ^AM = 6.5 X lO-SMoyr-i) 


0.0915 


1400 


To,p{M ^ 1.625 X lO-^Afoyr-i) 


0.0802 


1402 


ro^^(A/ = 3.25 X 


0.0801 


1453 


Ta,p{M = 6.5 X lO-^Mgyr-i) 


0.0971 


1430 


Eisner et al. (2009) 






RY Tau 


0.16 


1750 


DG Tau 


0.18 


1260 


RW Am- 


0.14 


1330 


AS 205A 


0.14 


1850 



mass star, the sublimation radius (equal to the wall location) is Rsub = 
0.079 AU . The height of the wall is ta ken as 5 times the pressure scale height 
([Dullemond. Dominik k, Natta 1120011) . Note that the parameters required to 
define a T^jix wall is Tsub (or Rsub) and a height (see Figure [T]). 

The SED of this model is presented in Figure [5J In order to do a fair 
comparison, wc add to the wall and star spectra, a model of a disk with an 
inner radius equal to the location of t he wall and an o u ter ra dius equal to 
IQQAU. The disk model is done using iD'Alessio et all |l998l) . The values 
of Tsufa and Rsub for this model are presented in Table [TJ Note that the wall 
spectrum does not show the lO/im silicate spectrum. Models changing the 
disk inclination and the mass accretion rate arc presented in § 13.21 



3.2. To,p WALL 

The shape of the T q.q wall is given by the fact that Tsub depends on density 



([Pollack et al. Ill994f ). The densest parts close to the midplane have a Tsub 
larger than the Tsub in the upper layers of the disk, thus, the former is closer 
to the star than the latte r. The density vertica l profile is taken from 2D 



axisymmetric disk models (jD'Alessio et al. II1998I ) 



Figure |3] presents a plot of Tsub vs Rwaii for this case, which shows a 
decreasing temperature with the radius, as one expects from a disk with a 
decreasing density with an increasing radius. The shape of the wall is given 
in Figure [1] either for Tg.p or To./ia: walls for three values of M. The value 
M = 3.25 X lQ~^MQyr~^ corresponds to the fiducial model; for comparison 
two models with half and twice this value are presented. Note that the T^jix 
wall location is given at the position where the temperature is equal to 
The Tq^p walls are located consistently outwards of the Ti^jix walls. 
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Fig. 2. SEDs for the T^ji^ model (thick hne), and for the To,p wall shaped by the 
dependence of the sublimation temperature on density (thin line) for the fiducial 
model. The dashed lines represent the wall spectrum, the pointed lines the disk 
SED and the solid lines the total emission. The disk inclination is cosi = 0.5. 

Figure [2] also presents the spectrum of the fiducial Tb,p wall model. The 
model includes the outer disk, starting in this case at the outer radius of the 
Tg^p wall. The T^jix wall emission is higher than the To.p for A < 8/xm. Also, 
the emission of the disk associated with the Tqjix wall is higher because in this 
case the outer disk start at a smaller radius. Putting these facts together, the 
SED of a system with a Tqjix wall is noticeably higher than the model with 
a To_p wall. This occurs mainly between 3 and 8/im, resulting in differences 
between the IRAC colors, as it is described in § |4l 

Neither the mass accretion rate or the inclination are parameters not well 
defined or even not defined at all for real systems. In order to see models 
with different values for these parameters and to be sure that the models 
do not overlap, we present another set of models. Figure |3] presents for the 
fiducial model, the behavior of the spectrum with inclination, either for T^jix 
or To_p models. It is important to note that for each type of wall, the emission 
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Fig. 3. The sublimation temperature Tsub along the To.p wall surface is shown 
with a solid line for M — 1.625 x 10^^ A'lQyr^^ , with a pointed line for M — 
3.25 X 10~®MQt/r"\ and with a dashed line for M = 6.5 x lO^^Mgt/r^^ 



decreases as the inclination increases. Note that for a higher inclination, the 
surface of the disk in the sky plane decreases, which naively means a lower 
emission. Besides, note that at A < 8/xm, a T^jix wall model with cosz = 0.5 
(60°) emits more than a To^p wall with cosi — 0.7 (45°). It is important to 
point out that the shape of the SED for these two models is different, thus, in 
principle a change in inclination is not able to match models with a T^ p and a 
To, fix wall, thus, in this way a modeler should be able to distinguish between 
these two scenarios, of course depending on the resolution and precision of 
the spectrum. 

In Figure \E\ plots for the SED in terms of M are presented. It is expected 
that either for the T^jix or To.p wall models, increasing M means a larger 
flux, and this do not mean just to move the SED by a constant amount, but 
a change in the shape of the curve. Note that the emission for a T^jix wall 
model with M = 3.25 x lO^^MQ-yr^^ is larger than the emission from a Tp.p 
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Fig. 4. The spectrum of a Tgji^ model (thick lines), and the To,p model (thin lines). 
The inclination taken corresponds to cosi — 0.3 (solid lines), cosz = 0.5 (pointed 
lines), and cosi = 0.7 (dashed lines). 

wall model with M = 6.5 x 10~^Moyr~^ for A < 7/im. This is an example 
of the differences between To,p and Tojix wall models, which correspond to 
changes in the IRAC colors as can be seen in § 21 

4. COLOR INDICES 



In I Allen et al. I (|2004[ ) there is a comparison of IRAC colors between mod- 
els with Tgjix walls and a sample of young stellar objects. Here, we note a 
difference in the SED between the models with Tojix and To.p walls in the 
wavelength range defined by 3 and 8/im. This means that the IRAC colors 
will show variations between the models. The emission of the Tojix wall that 
they use is a blackbody at the typical sublimation temperature of 1400K. Note 
that the emission from our model for the T^jix wall comes from the atmo- 
sphere, with the surface at 1400K, but the inner parts a little bit cooler. This 
means that ours and theirs colors for models with the same parameters should 
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10 
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Fig. 5. The spectrum of a Toji^: model (thick lines), and the To,p model (thin 
lines) in terms of mass accretion rate. The mass accretion rate taken corresponds 
to M = 1.625 X 10"*MQi/r-^ (solid line), M = 3.25 x W'^Ad^yr'^ (pointed line), 
and M = 6.5 X 10~^MQyr~'^ (dashed line). 



not b e the same. The [3.6]-[4.5] color range coincides between lAUen et al. 



(|2004l ) and this work. The [5.8]- [8.0] color value associa ted to our models is 



consistently larger than the values in I Allen et al. I (|2004l ). Note that even for 



the T^jix walls, the emission comes from material in a range of temperatures, 
in parti cular at temperatur es lower than 1400K, which is the temperature 
taken in lAUen et al. I (|2004[ ) . This should change in particular the [ 5. 8] -[8.0 ] 
color. The zero point magnitudes are taken as in Hartmann et al. ( 2005bl ). 



consistent wi th a Vega-based I RAQ magnitude system 



Recently, iMcClure et al. 



(|2010[ ) estimates spectral indices in the range 
between 6 and 31/im for disk models, using Tq ji^ walls at 1400K. As noted 
here, a change in the emission due to a To^p wall instead of a Tojix wall 
wi ll modify the spectral i ndices, however, due to the dispersion of the values 



McClure et all (l2010l) . including a To.p wall will not change their main 



14 



NAGEL ET AL. 



conclusions. On the other hand, when one tries to fit specific objects, taking 
into account a Tg.p wall should be an unavoidable modeling requirement. 

In Figure [6l the IRAC color-color diagram is presented for Toji^ and 
To,p wall models. Models with M = (1.625,3.25,6.5) x lO'^M^yr-^ and 
cosi — (0.3,0.5,0.7) are presented. The cosi ~ 0.3 Toji^ wall models are 
clearly located around [3.6] - [4.5] ~ 0.4 and [5.8] - [8.0] ~ 0.65. The less 
inclined models moves up and right in this plot, either for the Tgji^ or Tg^p 
walls models. The To_p wall models compared with respect to their T^jix 
counterparts, decrease the [3.6] — [4.5] color and increase the [5.8] — [8.0] color. 
The trend of the colors in the plot changing M and cos i, allow us to confidently 
say that there is no overlapping of the To.p and To./m wall models. Thus, 
one cannot confuse a T^jix and a Tg.p wall model when changing these two 
parameters. 

In the last paragraph, we stress that the various models are located in 
specific regions of the color-color diagram. Specifically, looking at the models 
in Figure [6l there is no degeneracy with the inclination. At first sight, this 
looks as a surprising result, because it is usually assumed that the inclination 
can work as a tuning parameter for the emitted flux. Notwithstanding, it 
is highly probable that we can tune one of the IRAC magnitudes, but it 
is difficult to think that one can tunc the four at the same time. In favor 
of this, this is not just a geometrical problem, because the fiux comes from a 
temperature distribution that depends on the position in the wall atmosphere, 
different between the models. 

In order to compare the modeled color indices with real s ystems, in Figure 



the bla ck points correspond to the Class II objects of the iHartmann et al 



3 



(j2005b[ )'s sample of Taurus. A direct comparison point out that the objects are 
consistent with high inclination disks. However, a definitive answer requires a 
modeling of the systems one by one, which is beyond the scope of this paper. 

Summarizing, the results presented in Figure [S] lead us to the conclusion 
that the models are not degenerate for the parameters used. In particular, if 
one can observationally fit the inclination (i.e. using interferometric images), 
a comparison of the observed and modeled colors will allow to choose between 
the models. 

5. SUMMARY AND CONCLUSIONS 
In this paper, for simplicity the heating from inner layers is not taken into 



account. This is responsible to move outwards the rim ([Kama. Min fc Dominik 



20091), and also to change the shape of the wall. Here, we describe the wall as a 



stationary well defined surface, in spite of the results of lKama. Min fc Dominik 



(|2009l ). They developed a code that find the sublimation wall with a detailed 
description of the sublimation process: the wall (dust-gas boundary) is a non- 
stationary diffuse region. Our reason for this assumption is that the calcula- 
tion of the emission for a time dependent diffuse region is a very complex issue, 
which depends on a lot of unknowns. For the lack of this knowledge, to pursue 
further in the improvement of the model is not worthed at this moment. Our 
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Fig. 6. IRAC colors for models including star, wall and disk. The models with a 
Tojix wall or a To^p wall are presented with open and filled symbols, respectively. 
Models for the following mass accretion rates are shown: M = 1.625 x 10~^ M^yr'^ 
(circles), M = 3.25 x lO'^Mgyr"^ (triangles), and M = 6.5 x lO'^Moj/r"^ 
(squares). The smaller size symbols correspond to an inclination given by cos i = 0.3, 
the medium size correspond to cos i — 0.5, and the larger size to cos i = 0.7. 
T he black points correspon d to the sample of Class II pre-main-sequence objects 
in lHartmann et al. I (|2005bl ). 



goal is to detect the differences on the SED, in particular for changes on the 
vertical geometry assumption, which we think are not conspicuously modified 
due to these assumptions. 

There is a difference between the Tojix and To^p walls SEDs in the wave- 
length range between A = 3 and 8/.jm. Due to this, the near-infrared colors 
calculated with a disk plus a T^ ji^ wall (commonly used) and a disk plus a 



T^ p wall differ. Howeve r, when analysing sets of spectra (I Allen et al. 112004 : 



Hartmann et al. Il2005bl ) using the IRAC color-color diagram, this difference 



do not change the conclusions previously presented, due to the dispersion of 
models and observations in this plot. For the modeling of real objects spectra, 
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the decision of which wall model to take, should be done carefully. 
The main conclusions are summarized next. 

• A To fix wall is closer than a Tq ^ wall. R^aii changes between T^jix or 
To.p wall models, at most by 10%. 

• The emission of a To,fix wall is larger than the emission of a To,p wall. 

• For each type of wall, the emission decreases as the inclination increases. 

• A change in inclination is not able to match models with a To.p a-nd a 
Tojtx waU. 

• The Tq jix wall do not show the 10/im silicates band (see figure [2) • 

• The disk for models with both types of walls shows the silicates band. 
For this reason, in the spectrum of the star-wall-disk system is always 
present this feature. 

• In IRAC color-color diagram, the Tqjix wall models are located in a 
region different with respect to the To.p wall models. Between the To,p 
and To, fix wall models, the [3.6]- [4.5] and [5.8]- [8.0] colors changes in 
around 20% and 10%, respectively. 
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